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l^tl-Oroup,  HultioHef lector  Pile  Theory 

1«  Introduction.;  The  differential  eq^^Ltlons  Involred  in 
the  taaltl»-eroup»  itfultl^refleetor  pile  theory  studied  In  tb^s  paper 
are  hejs^gencous  partial  dlff  erential  etpxatious  of  the  second  order.. 
Asfio elated  with  the  pile  or  any  of  the  refleetors  ie  a  ^sten  of 
at  many  e^juatlons  afi  there  are  groups<.  Thus,  if  there  are  a  iproups 
and^  reflectors  there  is  a  system  of  n  second  order  equations  la 
each  of  the^>4'  1  regions,  ths  solutions  of  the  systeas 

satisfying  appropriate  ooatinuity  conditions  at  the  boundaries. 

There  are  three  types  of  pilss  eonsidersd:  the  infinite  slab» 
the  infinite  dyllndcro  and  the  sphere.  Moreover,  syonetrical 
characteristics  of  the  problens  are  such  that  the  partial  diffemtial 
equations  reduce  to  ordinary  differential  equations^  Although  the 
equations  are  linear «  the  coefficients  are  constant  only  in  the  case 
of  the  slab  pile^  otherwise,  they  are  fTinctions  of  the  independent 
variable. 

Ordlnsprlly,  the  problem  to  be  solved  is  to  fixui  the  neutron 
density  at  any  point  within  the  region  bounded  by  the  outermost 
reflector  and  to  find  the  critical  size  of  the  pile,  after  the 
diaensions  of  the  reflectors  have  been  assigned^  In  all  three  cases 
mentioned  the  general  solution  of  the  set  within  each  region  ie 
easily  written  down,  even  where  the  coefficients  are  variable^ 
Thereafter,  a  difficulty  arises  in  that  the  location  of  the 
boundaries  is  not  assigned  but  is  dependent  upon  the  else  of  the 


_-lf= 

pll9.,  ISbne,  A  procedure  for  the  eolution  of  the  prohlea  ie  to 
soItc  the  eqaatto&o  aequeutlally  from  the  pile  through  all  of 
the  reflectors  in  terns  of  a  fixed  hut  unspecified  critical  else 
(radius  or  half »width)  until  the  outer  boundary  Is  reached, 

ThoBt  the  outer  boundary  oonditlone  tog^ether  with  re^ireneets 
of  synnetry  within  the  pile  proTlde  sufficlaxt  speolfications 
for  the  determination  of  critical  pile  sine  and  the  coaplete 
solution  of  the  probleo. 

In  a  prerloue  report,  Mon  P-^202,  the  two-group  theory  for 
the  ease  of  nultiple  reflectors,  without  oultlplication  in  the 
refleotore,  and  an  Infinite  slab  pile  is  dleeussed.  In  the 
deteralnatlon  of  critical  else  it  is  found  that  the  contribution 
of  each  reflector  to  the  solution  of  the  problea  ie  independent 
of  pile  else.  The  equation  which  finally  yields  the  critical 
size  Involves  a  single  fourth  order  determinant,  whoea  first  two 
columns  only  involre  the  pile  else,  and  idu>se  last  two  columns 
are  obtainable  as  a  product  of  0=1  fourth  order  matrices  and 
one  Uk2  matrix  (four  rows  and  two  columns),  all  independent  of 
pile  slxec 

In  this  paper  the  problem  is  discussed  from  a  more  general 
point  of  view  which  includes  the  consideration  of  possible 
staltiplication  in  the  reflectors  and  throe  different  pile  shapee, 
?lret  of  all,  the  system  of  n  second  order  differential 
equations  is  replaced  by  a  system  of  2n  first  order  differential 
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aquations.  This  procedure  eabraees  all  of  the  adTanta^es  of 
the  first  etndjr  while  further  providing  better  general 
perepootlve  on  the  problen  and  sloplifylng  soaewhat  certain 
oos^tttational  aspects  of  the  problea.  The  slaplieltj  first 
achieTsd  In  the  Infinite  slab  case  is  preserved  in  the  new 
procedure,  but  does  not  persist  In  the  spherical  and 
cylindrical  cases  for  either  siethod.  However,  In  all  eases 
it  is  possible  to  reduce  the  order  of  the  fundaaental 
deteminant  for  the  detersinatlon  of  critical  pile  else  to 
exactly  the  noaber  of  groups  involved  in  the  problea. 

In  ^  2  a  epecial  exaaple  is  provided  for  the  ease  of  two 
groups  and  an  arbltrsury  number  of  reflectors  in  an  infinite 
slab  pile,  with  aoltlpli cation  In  the  reflectors.  In  |  3  the 
case  of  the  Infinite  slab  pile  is  treated  with  full  generality, 
Oenersl  treatments  of  the  spberloal  and  oylindrioal  eases  are 
given  in  ^  h  and  ^  5  respectively,  finally,  in  ^  6,  a  study  Is 
made  of  the  limiting  situation  in  the  case  of  the  Infinite  slab 
pile  and  the  spherical  pile  to  provide  on  n  »  group  theory  with 
a  reflector  of  variable  densltyo 

2”  The  tw^group .  anlti^ef lector  problem  for  the  infinite 
slab  pile. 

2.1.  the  differential  e goat ions  j^d  boundary  eondlttone. 
Let  y,  with  an  appropriate  subscript  to  represent  a  particular 
group,  represent  the  flux  of  the  neutrons  of  this  group  in  any  of 
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the  1  regions  oonelstlog  of  the  pile  and  the  reflectors. 
Let  7^  represent  the  flux  of  thermal  neutrons  and  the  flux 
of  non^thermal  neutrons.  The  equations  in  an7  region  have 
the  fora 


(2.1.1) 


dx^ 

2 

»  72 
dx® 


where  the  constant  k  is  sero  In  a  reflector  without  ■ultlplloatlon. 

The  boundary  conditions  prescribe  the  continuity  of  y|^  and 
dsT 

1  ,  for  1  =  1,  2,  at  the  Interaedlate  boundaries  and  the 

dx 

vanishing  of  y^  at  the  outer  boundary. 

2.2.  Heductlon  to  a  system  of  first  order  differential 
equations.  If  in  (2.1.1)  the  substitutions 


^3 


y 


4 


y 


■j 

2 


are  aade,  the  system  of  first  order  equations  which  follow  is 
i^bialned^ 


Soltt^loa  of  tho  proMea  for  tfao  case  of  ao 
aultlplleatloB  la  the  reflectors.  Solutioas  of  the  fora 

s  e'-'  *  •  1  a  1,  2,  3,  4,  are  aoe  sou^^t.  Substitution 
la  (2.2.1)  yields  the  equations 


.8. 


-1 

>2  *4  • 

44AS  “  k  \''^1  ^2'‘'2  h  * 


which  hare  a  noc^trlrlal  Bolntion  In  A^,  Agt  ^y  provided  thei^ 


<2«3.1) 


0 


\ 


/*- 


>.^-y~l 


\=.l 
A  2 


0 


0  • 


Thle  Is  the  characteristic  eqnation  of  the  s^stea  (2^2,1)., 
The  equation  (2. 3,1)  can  he  written  in  the  fora 


0 

0 


0 

0 


-Vl 
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aad  coB«eqa«nt]y 


or 

(2.3.2)  CKf  (T<|  -^2)  -  k  'X|  s  0  , 


In  tbo  pile  k  >  1  and  the  eqnation  (2.3.2)  In  /c.2  txae  one 
poeitive  and  one  negatlTe  root.  iRme,  the  four  roote  of  (2.3,2) 
Bay  he  written  In  the  fora  1:  1  'tA-2*  *'^*'**  /‘'g 

reel  and  poeltlre  nnahere.  In  thie  eeetlon  It  ie  aeeumed  that 
k;  a  0  in  the  reflectore,  and  the  eorreenonding  roote  of  (2.3.2) 
are  either  oaee  there  exiete  a  eet  of  four 

independent  eolntlone  of  the  eyeten  (2.2  1).  Let  each  eolation 
of  any  independent  eet  of  solutione  he  Bade  a  eoloan  of  a 


4  4  oatrix  T  ,  To  illuetratet  a  matrix  jgt  eorresponding  to 

the  »c  >  th  region  of  the  pile  ie  const raeted.  One  nses  a 

%  X 

solution  of  the  fora  y^  »  C  *  •  taking  ^  1.  in  the 

eqoations  (2.2.1)  to  get  the  first  ooluon  In  .  In  this 

..'Ki* 

case  it  is  fotind  that  Aj  -  Aj^  s  0  *  euhstitution  y^  *  A^C. 

glTss  the  second  eoloan,  and  so  on.  The  final  resalt  for  ^  >  0  is 


« 

Itotil  further  notice  id  this  section  a  suhseript  shall 
henceforth  designate  a  parti solar  region  inetead  of  a  parti colar  gronp. 


0 


>  K  • 

1  1 


2  12 


>v^X-  (X 


‘1  “^2^* 

■Kgx 

-w' 

X^)6  ^ 

•K'^CK*  -K. 

which  Is  a  aon-*iiiealar  utriz,  Scr  talcing  appropriate  linear 
comhinatlooe  of  the  column*  of  thi*  matrix  one  may  replace  the 
ahoee  matrix  by  a  new  one 


i  *  (»)  ” 


ehX^x 


0 


■h  z 
2 


e^ehXgX 

X^-KgChXgX 
Xg  s^XgChXgX 


ohX^z 

0 

\  -v  ahK  » 
^11  1 

0 


ohKgX 

Ch'KgX 


Xg  ej^Xgsh  ^ 


•K®) 


lAere 
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for  later  reference  the  matrix  T^,  corresponding  to  the  pile, 
ie  computed  to  be 

^  (*)  » 


cos/<-^x 

chAgX 

•i“Aj* 

shAgX  \ 

r^^eos  /A.^x 

2Ai*^/‘'i* 

r  X  jif  sh  X 

2  2^2  ^2 

'2^  2'^2**A2*  / 

tdiere 

/ 

\  \  '*2  2  2 

ri®\l>2T<2  • 

V  -1  -2  2  2 

'2“\^2  *^2  <”’^1  "Aa)  » 

One  heart  in  mind  that  the  parameter*-;  and  >Mdilch  occur  here  are 
those  characterising  the  material  of  the  pile  and  are  distinct 
from  those  oeeurrlng  in  the  matrix  Tx  written  aboT*, 

In  general,  then,  the  matrix  T(x)  for  each  region  Is  a 
non-singular  matrix  of  functions  of  x  satisfying  the  matrix 
^fivation 

T*  «  W. 

Since  N  Is  a  constant  siatrlx  the  related  system  of  e({uatlons  haTe 
constant  coefficients,  it  is  then  clear  that  if  in  the  matrix 


T(i)t  the  indepeEdent  variable  x  is  averywhers  replaced  by 
*  ..  Xgt  the  resulting  matrix  J{t  ^  ie  also  a  matrix  of 

indepesdeut  solutions  if  designates  a»y  point  vlthia  or  on 
a  boundary  of  the  region  (pile  or  reflector)  over  which  the 
differential  equations  arc  defined.,  Indsedc  the  rsplaceraent 
of  £  by  X  ..  x^  is  serely  a  traneforaation  of  the  independent 
variable  which  has  no  effect  upon  ths  form  of  the  equsttoaso 
Moreover,  every  solution  of  the  eqiistlone  is  expreeslble  In 
the  fo.rs  y  Yc.  yhero  c  is  a  colu.'jin  vector  of  four  elenents.- 
Oa®  now  takes  the  origin  as  the  center  of  the  pile  t^nd 
X  and  x^  so  the  bounds  of  the®^  -  th  region..  Accordingly* 
the  aiatirlx  be  ^aicea  as  a  fwaction  of  x  «  x  e  le^ 

of  th®  difitaxic©  from  the  Inner  boundasry  of  th©  re^on  to 
arbitrary  point  of  the  regiosio 

The  conditions  of  continuity  at  the  pile  and  reflector 
boundaries  require  that 

• 

or,  if  c  x^-  X  , 

y^  (t^  )  ^  y^^^  (D)  o 

notice  that  in  the  last  equation  neither  x^  nor  x  ,  appears* 
but  only  the  tbicknees  to^.  of  the  *<.  -  th  region. 

Suppose  now  that  a  solution  y^  were  given  for  the  pile. 


This  could  he  expressed  In  the  fora 


y©  s  ^0  ®o 

idiere  Tg  is  the  particular  siatrix  of  fundamental  solutlohe  whldi 
has  been  chosen  vlthln  the  region  of  the  plle<  low.  the  solution 
Tf\  Ch^  in  the  first  reflector  aust  satisfjr  the  relation 

irtiere  a  s  tg  is  the  balf>thidcness  of  the  pile. 

Hence. 

and.  since  all  of  the  '■  hare  been  chosen  non->singDlar.  one 
obtains 

c^  a  T“^(C)Tg(a)e^  , 

Again,  at  the  next  boundary. 

ygCo)  s  yjCbj)  . 


12(0)02  a  Tj(tj)0j  , 


or 


-Ill- 


or 

If  thio  proeodnre  !■  eontinaed  oae  obtalu 

(2.3.4)  o^  “  °  • 

Thuti  if  iho  Toetor  (or  tho  tolution  vlthin  tho  pilo)  and 
tho  balf-thldcaost  a  of  tho  plla  wera  known,  feha  aolutloa  within 
all  of  tho  refleetora  eould  ba  obtalnad. 

Iha  aitnation  laat  eonaidarad  doaa  not  ganerallr'  ariaa,  and 
both  and  a  mat  ba  obtalnad  froa  othar  eonditlona.  tha 
phjraleal  aitnation  radpiraa  that  tha  aolntion  ba  arnmatrie  within 
tba  pila.  deoordinfily,  tha  laat  two  eoaponanta  of  tha  wactor  Jq 
mat  wanlah  at  tha  outar  bonndaxyi  that  ia,  tha  firat  two  eoi^wnanta 
of  tha  wactor  mat  waniah  at  tha  outar  boundary.  To  azpraaa 
thaao  raduiranenta  in  aynbolie  fom  lat  and  Ty  ba  partitionad 
into  2x4  natrleaai 


and  than  raqulra  that 

(2.3.5)  V^®>®0  «  °  • 

(2.3.6)  »^j(t^  )co  s  0  . 

If  tha  matrix  1»  wrlttm  in  tha  fom  (2,3.3)  tha  radoiramant  (2.3.5) 


=.15=^ 

l8  met  tgr  oqkiog  the  laat  ttK)  components  of  the  rector  eijual 
to  sero.  The  end  result  is  the  same  if  one  drops  out  the  Inst 
two  columns  of  thus  mahing  It  a  4  x  S  matrix  cooroosed  of  the 
first  two  columns  of  (2,3>3)»  end  it  will  ho  asstindd  that  this 
has  been  done  in  what  follows^ 

Bow,  In  (2,  3<4),  fiet  oc  »  J  and  maltlpl;  on  the  laft  by 
j^(tj  )„  Then  It  follows  from  (2,3»6)  that 

■^1  1 

The  matrix  I  (0)y  (t  )l“  ,(0),  =  ,  Is  a  4  *  2  matrix 

V  V  V  --'X  ® 

since  T^(a)  le  a  4  x  2  matrix,  and  ^(to  )  is  a  2  x  4  matrlXc 
Then  the  product  is  a  2  x  2  matrix.,  The  equations  are  consistent 
if  and  0.UI7  if  the  determinant  of  the  coefficients  ranlshes; 

(2^3.7)  ,  I^(a)\  »0  , 

The  determinant  is  of  order  two  and  the  solutloh  of  the  equation 
gives  the  critical  pile  else  a  . 

The  computation  is  simplified  in  a  numerical  nrohlea  if 
Y.,  (x  2  , )  is  replaced  by  T  j  (x  -  x  1  >„  This  has  the  effect 

of  sjnfelng  th<9  two  comronents  of  the  vector  c  ^  eijvusJL  to  serOfl 

find  condition  (2o3<^o)  becomes 

I-V>i(0)c's0  » 


equation  (2,3-7)  is  oow  written 


The  first  two  rows  only  occur  In  the  dBternlnant  of  (2«37)- 
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It  should  he  emohaeiaed  that  in  the  preceding  formulas  the 
subscripts  on  the  Y's  and  t's  refer  to  a  particular  region  (pile 
or  reflector)  and  the  subscripts  elsewhere  designate  a  particular 
group.  When  these  formulas  are  used  for  actual  computation  a 
double  subscript  notation  (IJ)  should  be  used,  where  i  deeignatee 
the  region  and  J  deeigxiates  the  group- 

2  4-  I'urther  observations  on  the  solution.  It  it  possible 

that  the  computation  In  the  problem  of  the  preceding  section  can 

be  simplified  by  making  the  matrlr  T  ^  (0)  of  Initial  values  the 
identity  matrix.  It  Is  clear  that  if  I  ,  (x  »  *  1»  W 

singular  matrix  of  solutions  then  so  also  is  the  matrix 

V  X  )2T,^(x«x  )Y.  (0)  , 

and  i8or«over 

t  ,  (0)  3  I 

Then,  the  fundamental  e<luation  (2.3'>7)  becomes 

1  .)  “  **02^*^  1  ”  0 

whoro  )  conEiBtg  of  the  first  tvo  rows  of  (t^  )  and 

consists  of  the  flr&t  two  columns  of  Wq(&)o 
In  the  preceding  section  the  formula 

(2-4.1)  po^ 

wa*  obtained  without  Imposing  any  restriction  on  the  matrix  of 
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Mlutlone  beyond  Ite  noxweloeal^rlty.  It  la  now  obeerTed 
that  (2.4.1)  can  be  written  In  the  eqolTalent  fora 

Cb  a  To^(n)Tl(C)Ti^(ti)  .  .  ,  Tg(0)ey  s 

Aeeonpanyiac  oonditlona  required  that 

(^2  •  3  0  * 

(Ig  •  s  0  * 

where  O2  and  Ig  are  2x2  aero  .  and  Identity  .  natrleee. 
reapeetlrely.  fly  the  special  choice  of  Tq  whereby  a  partitioning 
of  Tq  (0)  has  the  fora 


the  condition  on  the  rector  alone  la  equlralent  to  the 
requirement  that  the  last  two  components  of  Cq  be  aero.  The 
aubaeqaent  elimination  of  the  two  rectors  Cg  and  wae  found  to 
yield  a  determinant  of  osrder  two  la  the  final  equation  of  the 
problem.  A  eltnllar  special  dliolee  of  Tg  la  useful  \dien  the  outer 
reflector  Is  Infinite  In  the  x  .  direction.  To  this  end  one 
ehooses  the  first  two  columns  of  Ty  to  be  a  set  of  solutions  of 
the  differential  equations  %dileb  ranlsh  at  Infinity.  Fence,  the 
last  two  components  of  are  both  aero  and  there  remains  the 


determlnantal  equation 

(2<>4<2)  °  ’  *  -  0  • 

eliere  le  the  lover  half  of  the  aatrlz  Tj^(O)  and 

th*  left-hand  half  of  the  natrix  T^Ca). 

A  procednre  idilch  elopllfiee  the  eoaputatlon  of  ^ngCO) 
In  equation  (8.4.2)  hae  bem  ensgeeted  by  B,  Splnrad.  A 
epeoial  oholee  of  the  laet  tvo  eoluBute  in  is  aade  eo  that 
the  taek  of  coniDUtln*  Tj^(O)  ie  made  ae  elnple  as  poeelble. 
One  writee  fpCx)  1&  the  fora 
T«(x)  . 


.Xt* 


X.'K  j^e 


.X2* 


X  2®^  *K  2^ 


•h  K 1 X  1 

0  5^. 


chX,x 


i.(di  k-x  -  ch 

•l 


so  that  T_(0)  tehee  the  fora 


A  0 


B  I 


the  inrerse  of  which  le  readily  computed.  It  ie  the  lower  half 
of  thie  inveree  whi^  ie  need  in  (2.46  2). 


20: 


In  case  the  pile  is  cot  eymmetric  there  Is  no  reOnirenent 
at  symmetry  at  the  center.  Instead.  If  0  and  5  df^eignate  the 
two  outer  reflectors.  It  le  required  that 


(Ig,  0g)Io(0)c<,  »  0  . 

(Ig*  ^  0  • 


proTldod  that  both  outer  reflectors  are  flnltSo  The  aneclal  c}-oice 
of  both  rjad  Yg  results  in  Y^(0)  and  Tjj(^)  having  the  partitioned 
fora 


0 

A 


This  means  that  the  first  two  columns  of  each  Tanish  in  their  first 
two  components  at  the  outer  boundaries.  In  c^jse  elt:,er  reflector 
is  infinite  in  the  x  ^  direction  the  requirements  on  T^CO)  or 
Tjj(tg)  are  meaningless*  but  in  this  case  the  first  two  columns  of 
the  T  matrix  are  chosen  to  vanish  at  infinity ^  In  either  case* 
the  last  two  components  of  c^j  and  are  seroo  Hote  that  if  the 
sero  «  reflector  is  infinite*  the  origin  must  be  moved  to  Its  Inner 
boundary.  Now*  if  and  denote  the  two-vectors  of  non->sero 
components  there  results  either  of  the  equivalent  sets  of  2  x  2 
equations  of  the  form 


^ol  “  ^  ^  ^ 

from  iidilch  one  obtains  either  of  the  eqaiv^deut  determin  rji t  ,l 
equations 
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e  0  •  H  =0  . 


If  it  is  asauned  for  the  e^e  of  eiaplioity  that  all  of  the 
reflectore  are  finite,  note  that  it  aakee  little  differenoe 
vhich  of  the  B  redlona  ie  left  with  ite  thldcnees  iindeterBined;, 
Che  \inknown  thidcnesa  t^  ie  involved  in  only  a  elngle  Matrix  of 
the  product  or  E,  This  a.-.trix  is  Broltiplied  on  the  left  >y  & 
constant  matrix  of  dlmensione  2  x  U  amd  on  the  rl<dit  hy  one  of 
dimensions  4^2^ 

2,5«  Kultlplicetioo  la  the  reflectors^  In  this  section  the 
case  of  anltipli cation  In  the  reflectors  ie  considered.  If  k  s  0 
in  the  characteristic  equation  (2,3^2)  there  results  the  foreeoine 

p 

situation  In  which  there  are  two  positive  v%lue8  of/c  «  tf 

p 

0<.k41  the  yur*  B  •:re  etill  positive  (ind  the  matrix  is  as 
heforsv  If  k  s  I.  there  ie  one  positive/c^  and  one 
The  matrix  takes  the  form 


Y^  s 


/ 

(k2  +  "k|)x 

.(■k2  \ 

* 

X 

e 

\ 

Xl-Kf 

(1<f+'H|)x 

re 

.(x2+k|)x  1 

re 

XjKf 

0 

Xj. 

se 

.se 

^(■K2t1<|)x  / 

Be  1 

e 
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or 


\ 

Xi-Kf 

-TT 


■h(xf  +  T<|)* 


t  ah(  K  f 


s  cii(Kj  +  'v:|)x 


a  ch(Xj  +  'H|)x 


Vixf 


VgX 


a  sb(xj+x|)x 
a  8hCKj-»'^|)X 


vhere 


Xl(Kf  +  Ki)  . 


If  Jc  >  1,  T  (x)  take©  the  fora  of  Yn(x)  ae  la  (2,3,3), 


3’  polti^group^  i8aIti>:.reflector  problem  for  infinite 

slab  gile^ 

3.1.  general  forgmiation^  As  In  §  8,  let  y,  with  an 

appropriate  subscript  to  represent  a  particular  group*  represent  the 
flux  of  the  neutrons  of  this  group  In  any  of  theO-^  1  regions 
consisting  of  the  pile  and  the  reflectorn  Jhe  subscripts  are  chosen 
so  as  to  increase  with  the  mean  energy  of  the  group.,  that  Is,  y*|^ 
represents  the  flux  of  thermal  neutrons,  y^  that  of  the  slowest 
group  of  fast  neutrons^  and  that  of  the  fastest  group  (the  fission 
neutrons).  The  differential  eqaatlcna  in  any  region  have  the  form 


«23- 

=  0  ,  (  i  <  »  > 

» 

«b«r«  the  eonatant  k  le  aero  1&  a  reflector  without  naltlDlieatiou. 

Wie  boundar7  condltlone  preaerlhe  the  eontluult/  of  and  of 
%  at  the  Intermediate  boondariea  and  the  Tanlehlng  of  y^  at 

the  outer  hnondazy  or  honndarieSo 

3.2.  the  Infinite  alah  pile.  For  the  Infinite  eleh  pile 
the  differential  equations  are  (‘specially  simple  in  form  since 
^  becomes  merely  the  second  derlvatlYe  with  respect  to  X  and 
all  coeffielente  of  the  y^  and  their  derlvatlTes  are  constant. 

It  le  conrenlent  to  introduce  the  additional  rarlables 


(3-11) 


“  »  111  nna 


(J  s  ») 


Kiiich  ja‘Jc?e  it  poBcible  to  replace  the  syetesi  of  n  eeeond  order 
equatlone  by  the  cn  firet  order  equations 


‘'ll  yj  • 


(Jin)  . 


(3.2.1)  .  (l<.)  . 


4 


If  the  symbol  y  without  a  subscript  denotes  the  column  yector  of 
components  y^^,  y^^  -  ,  o  .  equation*  can  be  written  In 


the  matrir  form 


(3.2,2)  y“  g  1^  ^ 

wh«re  the  2n  »  2a  aatrlx  h  eaa  be  partitioned  in  the  fora 


Tb®  eharacterietic  equation 

)  fc  -  youl  \  s  0 


It  it  olaar  that  when  k  s  0  tha  roots  of  tba  eharaetarlstlc 
equation  are  To  iovestigata  the  eata  k>0  ecneldar  tba 

srapb  of  tba  left  mambar  of  (3»2, 3)  plotted  agalnet/*.  .  Tbit 
erottat  tbaA^  *■  axis  atK^  and  crosses  the  vertical  aria  at 

The  effect  of  tubtractiag  k  ,  .  ,  "V.  ^ 


t>IO 
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froffl  th«  fimction  Is  to  raise  the -  axle  hy  thie  aaoimt, 
Hsnce«  if  k  <.  1  the  real  roots  of  the  eijuatloa  in/*-®  are  all 
positive;  if  k  a  0  there  is  one  sero  root  and  all  other  real 
roots  are  positive;  if  k  > 1  there  Is  one  negative  root  and  all 
other  real  roots  are  positive*  Hovever,  one  cannot  be  assured 
in  general,  for  n>  2,  that  the  roots  In are  all  real  since, 
for  any  k  >  0,  by  making  the  difference  between  any  pair  of  'X 
eufflciently  small  one  may  introduce  a  pair  of  comulex  roots. 

Whatever  the  nature  of  the  characteristic  roots,  one  can 
always  employ  them  to  write  down  a  set  of  2n  Independent 
solutions  of  the  differential  equations.  If  each  of  these 
solutions  Is  made  a  column  of  a  2n  x  2n  matrix  T,  as  la  ^  2, 
then  T  is  a  non-singular  matrix  satisfying  the  matrix 
equation 


t’  s  Ksr  . 

^rom  this  point  on  the  discussion  for  the  case  of  tv» 
groups  given  in  |2  applies  without  significant  changes  for  the 
case  of  n  groups.  The  equation  for  the  determination  of  critical 
pile  Else  for  all  of  the  situations  considered  will  involve  a 
determinant  of  order  exactly  n  .  The  fundamental  determinantal 
equation  may  of  course  be  used  to  find  k  if  the  pile  and 
reflector  dimensions  are  fixed. 

multi- ref  lector  problem  for  the 


spherl cal  pile. 
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4<,  1 .  Tha  dlfferenti^il  etpiationa  unci  bonadary 
conditions »  If  the  epberic^il  pile  has  reflectors  in  the 
forn  of  spherical  shells,  sneb  as  to  provide  eosmlete 
sphericul  syaaetzy.  then 

r  'y  , 

dr" 

Then  the  substitution 


Vi  .  ryi 


^n  +  1 


.  9 

^i  ^ 


(i  i  »)  . 


in  the  differwtial  equations  (3-1-1)  yields  Sn  differential 
equations  in  Vi,  ^  i  of  exactly  the  sane  forn  as  those  in 
y^t  7j^^i  for  the  slab  pile-  However,  in  the  spherical  case, 
the  boundary  conditions  which  iAposn  continuity  upon  yi  and 
^1^1*  expressed  in  terms  of  Vi,  require  the  continuity 

of  Vi  and 

A  .  iii  _  Ali. 

‘  to  ■  I  ,2 

at  the  pile  and  reflector  boundaries.  Vote  in  passing  that 
idien  the  '>Si's  nve  equal  the  v^^i's  are  also  continuous  and 
the  procedure  is  in  all  respects  sidiilar  to  the  slab  ease. 

4-2.  Solution  for  the  eeneral  ease.  The  formation 
of  a  fundan(>ntal  set  of  solutions  v  constituting  a  non-singular 


matrix  7  does  not  differ  In  any  respect  from  the  formation  of 
such  -i  matrix  T  in  the  case  of  the  slab  pile*  These  solutions 
can  be  expressed  as  functions  of  r  or  of  P  m  r  «  r  «  uhere 
o<.  designates  the  x  «  th  region,  since  any  fundamental  set 
expressed  as  functions  of  r  remains  a  fundamental  set  uhen  r  is 
replaced  throu^out  by  In  this  procedure  it  must  be  noted 
that  one  obtain  solutions  y^  -  T^/r  which  become  infinite 
under  certain  cir corns tances^ 

It  is  suprK3sed  now  that  in  each  region o<  a  matrix  of 
solutions  Vx  has  been  found  and,  for  the  sake  of  conTenience, 
that  this  is  so  chosen  that 

Since  none  of  the  r^  ^s  known  urtil  the  critical  radius 
*^0  ®  found,  this  choice  inrolTes  expressing  the  solutions 

▼  as  functions  of  /o  ^  r  =.  r  ^  and  choosing  to  be  the 
identity  matrix  when  ^  0.. 

In  the  pile  (o<  s  0)  y  must  remain  finite  at  r  ■::i  0,  so  that 
out  of  the  2n  solutions  n  can  be  eliminated  Accordingly,  the 
required  solution  is  of  the  form 

’'o  “  ^01  ®0 

wh«re  is  an  n  »  Tector  of  eonetant#  «jd  ie  th®  2n  x  a 
matrix  of  t» a  (axr^reeeed  la  tsraa  of  trigonometric  and  hyperholl 
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tine*)  Uxoae  first  n  rows  waslsh  at  r  «  < 

T  ,  onst  ts  A, 

Oj. 

fhs  oontlAolty  re<ialreBsnt  lAToXyes 
vector 


(4.2.1) 


1 

0 


0 


0 

I 


0 


0 


>^2 

«»■»»» 

r^ 


0 


0 

0 


1 

0 

0 


0  0 

0  0 


0  0 


1. 

r 


0 


0 


1 

r 


0  0 


ThaSs  It  is  first  of  all  required  that 


h 


,  The  raBk  of 


the  eootlOAoas 
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providod  th<>»  s  1  a%  r  s  St 


Bote  that 


or 


I  ^2 ^^2^ ^2  “  ^3  * 

'  / 

vfaloh  is  usad  to  ellsdoats  is  the  last  t«o  aquations^  If  this 
proofldora  Is  continued  one  ohtalna  finally 


where  |  s  a+t^+tg  +  „  ,  „  vt^ 


This  mtriz  equation  expresses  the  2n  coaponents  of  In 
terns  of  the  n  conponents  of  and  the  xmknowii  critical  radltxs 
A.  An  additional  matrix  relation  is  needed  for  the  elimination 
of  e^,  and  this  is  found  in  the  requirement  that  the  solution 


shall  Tanisb  in  its  first  n  eoamonents  at  the  enter  houndar/o 


»32> 


Bov,  T  3  7^  .  nhence.  if  7^^  denotes  the  antrlx  ot  the 

first  n  rovs  of  7^  ,  the  condition  in  qaestion  takes  the  fora 

)®».’  s  0 

Aceordin^Ljr,  if  the  eq^ostion  (4,2o4)  is  smltiplied  on  the  left 
by  7^  there  results  the  eqnation 


The  matrix  coefficient  of  is  n  x  n  and  its  determnaot  anet 
Tanish,  that  is» 


Hot©  AgalsK  that  it  ha»  t^eea  presupposed  that  7^  (0)  ^  if 
this  xiosraalisatioai  is  not  made  the  matrices  To(°^(0)  mast  also  he 
included  lu  the  aatrlx  product Hot©  also  that  a  does  not  occur 
la  aay  of  the  7^  s  foT^?  0  hut  only  in  the  elements  of  the 
fcim  -vt^  ),^  multiplied  hy  a  difference  of/v^e  in 

coneecutire  regions.  Bote  finally  that  7^  (t^  )  is  of  the  sai”e 
fora  as  in  the  slab  ease, 

5-  SS  ^lti,=group.  BultUrafleetor  problem  for  the 
infinite  qsrltndrlcal  plleo 


5^1.  Tha  dlfforoatial  «qa»tloa»  ytd  boondary  oondltlone, 
It  the  cylindrical  pile  hae  reflector*  In  the  form  of  oarlladrical 
Shells  each  as  to  proride  axial  eymmetry,  then 

o  d  ^  1  dy 

^  7  -  p  r  dr 

dr^ 

The  differential  oqoatlon*  in  any  region  nay  ho  written  in  the 


form 


h  \^l  ^i-»  1^  1+  1*^1-^  1  '  ® 


(1  <  n) 


^'n  i-  (r  xfry  -K^  ® 

dr  dr  111  n  n  n 


Th©  ’bpiindAiy  ooAditiooG  r©^ird  the  continuity  of  JTj^  ©nd 

>V,  at  the  intermediate  hoondariee  and  the  Tanlsbing 

’■  dr 

of  yj^  at  the  outer  houndary  or  houndariee-, 

5-2.  Solution  for  the  general  case.  The  introduction 


of  the  variahles 


V 


(J  i  n) 


has  the  effect  of  replacing  the  eyeten  of  n  second  order 
Bipiatlons  (5.1.1)  hsr  the  2a  first  order  e<iaatlon« 


ri  .  \^1  _“i 

y^  =  r  y^^i 


(5.2  =  2)  y^^js  - 


(i  i  n). 


(i  <  »)» 


«34= 

If  the  ^nhol  y  without  a  Dubacripi  deootes  the  colusia  xrector 
of  coapoaeate  yg,  c  »  ,  ,  y^^,  thea  the  oquatloaa  (5o2o2) 
eaa  he  wlttan  la  the  matrix  form  y"  e  My,  where  M  le  the 
Boa^eonetaat  matrix  of  coeffleieate  la  the  eqoatione. 
Solutloae  of  the  aystea  (5.2„2)  of  the  fora 


<5.2.3) 


^1  ^  » 

are  now  Bought^  Suhetitution  of  (5.2^3)  lato  (5.2.2)  glvee 

(5.2.4) 


V-  n  a). 


(1  i  n). 

(1  <  a). 

These  equattone  have  a  non-trivial  soJutloa  in  provided  that 

-1 

I  1  ^  V 

(5  2o5} 


v«- 


yv 


^  0 


vhere 


0  .y;’- 


o  n  z 


<.'»  c 


0 


0 


\  ’ 
\o 


0 


-35= 


By  a  reduction  olmllar  to  that  used  in  |  3.2  the  ci'aracterletlc 
equation  t©ooffl«« 


(5^2.6)  n<f+A^) 


C  K^  §  feKJ 


Bote  that  the  root*  of  this  equation  Xu  /c^  ax’®  the  negatlree 
of  the  roots  of  (5.2,2) 

If  Jq  and  Jj^  are  replaced  hy  Tj,  and  respectively  In 

(5-2,5),  the  eyeteo  (5,2,4).  flJid  consequently  the  equation 

(5-,2-6),  reoalns  unchanged.  However,  If  the  comblnattone 

1q,  I,  or  are  used  the  effect  1»  the  expected  one  of 

2 

changing  the  si&is  of  the  roots  of  (5«2.,6}  in/<-  or  of 


replacing by  1/^  c 


The  notation  for  Bessel  fxinctione  need  In  this 
paper  is  concletent  with  that  fotind  in  Theory  of  Beaeel 
Functions  by  Q  B 


(5.2. 
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If  ttiG  roote  of  (3.2.6)  oro  all  jTiialt  thett  tho 

fonotlona  anfl  oorreqpond  to  tho  posit  Its  roots  and  ths 

fuoetions  1^  ocd  to  tho  nssatlrs  roots.  Thao,  in  this  eass 

at  loASt,  a  sot  of  2a  iadepoadent  solutioas  of  the  differoattal 
oipiations  can  to  foand.  4s  In  tho  preeedin£  ssetions*  if  ea^ 
of  thsBo  solntlone  is  nads  a  ooloam  of  a  2a  x  2ir  matrix  T, 
than  F  is  a  non-siagolar  matrix  satisfying  the  matrix  oqaatioo 

T  ®  S  M?  . 

In  tho  slat  and  sphorieal  eases  certain  material 
siatolifioations  were  found  hy  sKprsssing  tho  solntions  in  tho 
oc.  •  th  region  as  functions  of  a  coordinate  originating  at  tho 
innermost  hoondary  of  the  region.  This  siaplifieation  is  not 
possihle  in  the  (grlindrleal  ease  since  ths  matrix  H  is  no 
longer  a  constant  matrix,  /icoordingly,  the  determination  of 
critical  pile  siso  oron  in  tbs  simplest  cases  hoeoaes  an 
extremely  arduous  and  uninTiting  task.  The  eqoation  which 
fls^-ds  ths  critical  pile  sice  may  bo  written  as  the  n  x  n 
determlnantal  equation 

^  >)  ^  io  -l^^V  -=1^  .=1^%  »1^  \t)  *.2^°  ’  '»^Oi(»)\  SO, 

®  a  + 1^+  tg+  ,  „  a  +  t  ^  ,  7^^  is  a  &  X  n  matrix  of  tOs 
expressed  in  terms  of  the  functions  and  :i^,  which  are  finite 
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for  r  %  0.  and  olioro  ^  la  an  n  x  2n  natrixi  tha  vqppar  half 
of  tha  matrix  .  Hota  that  all  of  the  natrleaa  In  (3=2.7) 

Inrolra  tha  plla  atxa  a. 

mltl»groap  plla  theory  \flth  contlnnooely  raiylng 
paranetara  in  tha  rafloetor. 

6.1.  Tha  caaa  of  tha  Inflnlta  alab  plla.  In  thla 
dlscaaaion  It  la  nov  aaanaad  that  tha  X  <  a  and  'X  • a  of  tha 
Boltl-group  pile  thaoxy  aquatlona  ara  eontinuoualy  rarylng 
paramatara.  In  §  3  bha  systea  of  n  aaoond  order  dlff arantial 
agoationa  aaa  raplaoad  with  a  eyataa  of  2n  flrat  order  aqidatlona 
vith  the  aid  of  the  aubatltntlon 

idvara  the  'Xj'**  are  now  fnnetiona  of  .  It  follova  than  from 
(3.2.1)  that  tha  eyatan 

(  Xiyj  )®  —  J^l ^1  +  1 -»-l  ^1  1  ^  0  »  (1  <.  n), 

raplaeaa  tha  original  syatam  of  n  aeeond  order  diffarantlal 
agnationa. 

A  oonvani^t  point  of  depeurture  In  thla  aaetlon  la  tha 
eonaldaratlon  of  tha  aaqaanoa  of  aeta  of  aqoatlona 


tdisro  Xx  colottn  Tector  of  2ti  eotoponeats  ropretentini^ 

a  solution  of  ths  equatlone  in  th«  9(. »  th  reflector,  and  M  x 
is  a  constant  atatrix  eosposed  of  the  parameters  character! etie 
of  that  region.  Heretofore  no  relation  between  the  M" s  for 
distinct  regions  was  aesumed.  The  solutions  jr^  ,  bowerer, 
were  required  in  all  cases  to  satisfy  the  continuity  relations 

Xx  (»x  )  =  X  » 

If  y  is  expressed  in  terms  of  the  fundamental  matrices  Y  then 
the  continuity  relations  take  the  form 

(6.1.1)  T.,(s^)c.^s  i<*x  )c  x  +  l  ’ 

It  Is  supnosed  now  that  the  pile  and  reflectors  constitute 
a  single  region  with  continuously  vaspylng  parameters,  or 
parameters  haring  at  most  a  finite  number  of  discontinuities. 
Th«i  tbs  several  sets  of  equations  with  constant  eosffieients 
are  r^laced  by  a  single  set  of  eqnations 

T®  s  K(x)y 

with  varying  coefficients.  Since  the  case  in  which  ti  has 
components  with  disoontinoities  presents  no  essential 
eonplioations,  it  will  suffice  to  consider  only  the  continuous 
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caMo  Intarral  from  pile  toundaiy  to  reflector  bonodatr 
le  brolcen  op  iato  elabe  of  thidcnese  siad  an  appsroziiDatloix 
to  the  oln^e  eet  of  e<liiatloao  %rlth  Tatylag  ooeffiolents  le 
aohieved  trith  a  se^joeoee  of  eete  with  eooetant  ooeffieieatoo 
Thle  eeqoenee  of  eete  has  the  setae  fora  as  the  one  prerioosly 
coaeidered  with 

K»*WU).  >»  . 

that  le  to  say  with  M  of  the  fora  trfcen  by  M(z)  when  the 
eoapo&eats  are  eraloated  at  aa  arbitraxy  bot  fixed  point  within 
the  <<  »th  region^ 

The  fundaaental  mtrlcee  T^.  are  dioeen  as  In  ^  2^4  eo  that 

Y  =1^  ® 

Then,  from  It  follows  that 

0  oC  ^  f  ot  X,  )  ®  ot  => 

The  aatrix  T  ^  (x  x. )  aay  be  expanded  in  powers  of  A  x  to  glTo 

Tx  <*«  +Ax) 

s  T  .  (x  ,)  +  Tx. 


sC 


»^“l 


n 
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then 

It  vlll  suffice  In  this  dlseassion  to  keep  only  the  first 
two  terse  of  this  eeriee. 

By  definition  y  (x)  ®  7^  o  Tb»n  it  follows 

that  y.,  ®  s  Ic,^  c  e.,  „  Hence,  if  y(x) 

represents  the  solution  for  the  entire  reflector,  then  y(x) 
is  identified  tdth  yj^(x^)  ”  where  designates  the  inner 
houndsry  of  the  first  reflector  of  the  thi^^ese  A  x,  Horeover, 
the  entire  approxioation  to  the  exact  solution  y(x)  evaluated 
at  is  given  hy 

S  (I-^M*Ax)(I+M^^jAx)o^^^ 


(1+ A  x)(I+  A  x),  o  o(I-+M  Ax)c 

4  1 


vdiere  c  c^„  or 


r  l.*^  lyi 


(3 

1* 


■*■2 

In  the  passage  to  the  Unit  hold  fixed  the  point  x.,,  but 


allov  thd  RQB^er  of  su^diTifllons  to  increave  without  liodt 
irtilla  A  X  a|>proachee  *ero.  If  the  euhecript  «<,  is  suppreseed 
one  obtains  In  the  Unit 

y(x)a  j^l-ej  M(pd|+j  j  M(  ^  )M(  |»)d  t/4i 

+  1*  H(^  )M{  p)M(  |")dj*di«d|+  .  .  .  j  c  . 

If  the  matrix 

<6.1,2/  K(x) 1+  *  M(|  )df^-  [  'MCpiK  p)d|'d|  +  o  ,  o 

-  *0  *0  -^*0 

Is  Int.  ■  duoed.  the  solution  hecoaes 
(6.1,3)  y(x)  5.  K(x)c 

where  K(x)  is  a  matrix  whose  columns  constitute  a  fundamental 
set  of  solutions  and  whldi  reduces  to  the  Identity  at  x^. 

The  solution  (6,1,3)  can  be  obtained  by  a  well-fcnown 
classical  procedure,  first  of  all  it  is  noted  that  the  original 
equations 

y»  s  M(x)y  ,  yC*o)  “  ® 
are  equivalent  to  the  set  of  equations 
y(x)  a  C4-  f  )d_| 

Jxq 


0 
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Tba«.  bgr  repeated  substltutloae.  one  flnde  that 

7(x)sc+j  +  M(^07(p)d{«^  d^ 

s  j^l-v.  M(pdnc+-  I  f^«{  I  )M(  {  «)7(  iOdl'dI 


«  0  o  e 

This  procedure  glret  rise  once  again  to  the  matrix  1.  It  le 
knoun  that  this  series  converges  to  the  solution,  for  an 
arbltrar7  c»  under  ver7  general  eondltlonst.  The  general 
matrix  K(x)  ot  solutions,  or  the  particular  solution  7(x), 
could  be  approximated  b7  computing  a  sufficient  number  of 
terms  In  the  series  (6.,1«2),  The  repeated  Integrations  may, 
however,  become  ve«y  Laborious  or  totally  Is^sslble  to 
carry  out  la  closed  form.  Humeri cal  Integration  Involves  the 
breaking  up  of  the  range  of  integration  Into  sublnterrale  and 
the  Introduction  of  Interpolation  polynomials.  It  may  be 
simpler  and  less  laborious  to  apply  the  slab  procedure  for 
obtaining  the  approximation,  since  the  slab  procedure  ie 
no  more  complicated  tdien  the  slabs  are  unequal  la  thickneee, 
it  should  be  especially  advantageous  trtien  certain  parameters 
vary  at  hl^ly  non-uniform  rates,  permitting  greater  tbldcneeeee 
idjcre  the  r^tes  are  all  low,  and  requiring  smaller  subdivisions 
where  tome  are  high. 

It  is  worthy  of  particular  note  that  this  whole  procedure 
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It  fult*  indapaada&t  of  tho  fora  of  M  and  can  bo  appliod  to 
anjr  oot  whatsoorer  of  linear  eqoatlono. 

In  the  problea  of  detereiinlng  the  wltloal  tlse  of  a 
pile,  aeeuaed  to  hare  a  unifora  Interior  region  aithia  a 
reflector  of  eontinoouely  raiylng  paraaetere,  if  T  ie  the 
total  thldcnees  of  the  reflector  the  foregoing  procedure 
/ielde  a  relation  of  the  fora 

(6.1o3)  y(a+T)  ^  K(a^T)y<a)  . 

In  this  relation  the  a  appears  explicitly  in  y(a}  ^  o.  bat 
not  at  all  in  either  y(a-t-T)  or  K(avf)*  the  a-t-T  being 
aerely  a  place  label,  that  is,  y(at7)  and  E(at-T)  in  (6.1.3)  are 
functions  of  T  aionso  The  y(a-vT)  taJces  the  place  of 
Tj  (tj  )cj  in  §  2.  MoreoTer, 

y(a)  s  Ta(a)<^  s 

if  T^(0)  s  I.  Accordingly,  (^.>1.3}  can  be  written  in  the  fora 
y(a  +  T)  a  I(a-v-T)Tg(a)y{0)  , 

where  a  occurs  in  Tq  aloM^  The  first  n  eos^nents  of  y(a-«-T) 
are  sero  and  the  last  n  coaponents  of  y(0)  are  sero.  Thus  it 
follows  that  in  the  2n  z  2n  aatriz  Z(A+T)T^{ti)  the  ainor  n  z  a 
aatriz  in  the  upper  left-hand  comer  is  singular.  This  is 
obtained  by  aultiplying  the  upper  half  of  K  by  the  left  hand 


half  of 
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5,2,  Th«  oaoo  of  th>  ■pherical  pile.  The  dloeaooioQ 
in  tho  preeodlog  section  applieo  alaost  without  (diange  to  the 
ease  of  the  spherical  pile  with  a  reflector  with  oontiaaouslj 
▼aryiag  paraaeters,  fhere  is,  faoweTer,  a  significant  difference 
la  that  the  tranaforaation  employed  ahore  to  obtain  equations 
with  ooeffieients  that  are  constant  within  each  si^erical 
shell  at  the  same  time  proTides  solutions  which  are  discontinuous 
la  their  last  n  eoaponents  at  the  boimdaries.  The  substitution 
that  proTides  continuity  at  the  boundaries,  rlz,  ^ 

yields  a  set  of  equations  with  coefficients  idiich  are  functions 
of  r  erw  within  each  onifora  shell. 

The  adrantage  of  haring  equations  with  constant  coefficients 
is,  of  course,  that  their  solutions  can  be  written  down  readily 
in  terms  of  exponentlalso  or  tri^nometrio  and  hyperbolic  sines 
and  cosines.  Howerer.  once  solutions  are  giren  in  these  terms 
for  the  equations  with  constant  coefficients  (with  discontinuities 
in  the  solutions  at  the  boundaries)  a  simple  transformation 
prorides  the  corresponding  solutions  of  the  equations  with 
rariable  coefficients  (with  solutions  that  are  continuous  at 
the  boundaries).  A  further  readjustment  may  be  needed  to 
provide  a  matrix  of  solutions  reducing  to  1  at  the  inner  boundary* 
but  this  can  always  be  effected,  from  this  point,  the  procedure 
is  formally  the  same  as  for  the  slab  ease,  and  the  spherical 
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Bbell  appzozliaatlon  to  tho  ease  of  oontlnuously  varying 
paraaetere  can  be  approxloated  by  eolvlng  a  finite  aoaber 
of  eete  of  equations  for  unlfora  sbelle.. 

In  eomdaelon  It  is  noted  that  the  fora  of  the  solution 
in  the  oyllndrioal  ease  does  not  lend  itself  to  the  type  of 
analysis  need  in  the  preceding  eeetlone  for  the  problea  of 
eontlnnouely  vaxylng  paraaetersc  At  the  present  tlae  then 
the  iqrlindrioal  ease  with  oontinuously  varying  paraaetere 
in  the  reflector  oonstltutes  an  nnsolved  problea. 
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